MIMO OFDM system

ABSTRACT

A MIMO OFDM system includes a plurality of space-time encoders for encoding respective data blocks with independent space-time codes. The transformed data block signals are transmitted by a plurality of transmit antennas and received by a plurality of receive antennas. The received data is pre-whitened prior to maximum likelihood detection. In one embodiment, successive interference cancellation can be sued to improve system performance. Channel parameter estimation can be enhanced by weighting the channel impulse response estimates based upon a deviation from average.

CROSS REFERENCE TO RELATED APPLICATIONS

The present application is a continuation of U.S. patent application Ser. No. 12/621,335, filed on Nov. 18, 2009, now U.S. Pat. No. 8,121,022, issued Feb 21, 2012, which is a continuation of U.S. patent application Ser. No. 11/380,531 filed on Apr. 27, 2006, now U.S. Pat. No. 7,643,404, issued Jan. 5, 2010, which is a continuation of U.S. patent application Ser. No 09/791,523 filed on Feb. 23, 2001, now U.S. Pat. No. 7,068,628, issued Jun. 27, 2006, which claims benefit of U.S. Provisional Patent Application No. 60/206,182, filed May 22, 2000, all of which applications are incorporated herein by reference.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH

Not applicable

FIELD OF THE INVENTION

The present invention relates generally to communication systems and, more particularly, to Orthogonal Frequency Division Multiplexed (OFDM) communication systems.

BACKGROUND OF THE INVENTION

High data-rate wireless access is demanded by many applications. Traditionally, more bandwidth is required for higher data-rate transmission. However, due to spectra limitations, techniques for increasing bandwidth are often impractical and/or expensive. In one known system, multiple transmit and receive antennas are used to obtain transmit diversity or to form multiple-input multiple-output (MIMO) channels. Multiple transmit antennas have also been used to provide diversity in wireless systems. Transmit diversity can be based on linear transforms or by space-time coding. Space-time coding in particular is characterized by high code efficiency and can improve the efficiency and performance of Orthogonal Frequency Division Multiplexing (OFDM) systems. The system capacity can be further improved if multiples transmit and receive antennas are used to form MIMO channels. Compared with a signal-input single-OUTPUT (SISO) system with flat Rayleigh fading or narrowband channels, a MIMO system can improve the capacity by a factor of the minimum of the number of transmit and received antennas.

FIG. 1 shows a conventional OFDM system 10 including subsystems for transmission and reception of data. A coding subsystem 12 encodes binary data from a data source. The coded data is interleaved by an interleaving subsystem 14 and then mapped onto multi-amplitude multi-phase constellation symbols by a mapping subsystem 16. In one particular embodiment, the multi-amplitude multi-phase constellation symbols include quadrature phase shift keying (QPSK) symbols. Pilot signals can then inserted by a pilot insertion subsystem 18 to estimate the channel at the remote subscriber unit receivers. A serial-to-parallel conversion subsystem 20 converts the serial data stream to a parallel data stream that is provided to an inverse fast Fourier transform (IFFT) subsystem 22.

The transformed data is converted to serial data stream by a parallel-to-serial converter 24. Cyclic extension and windowing can be added by a subsystem 26 prior to digital-to-analog conversion by a DAC 28 and transmission by an antenna system 30. The receive portion 32 of the OFDM system includes similar corresponding components for extracting the data from the received OFDM signal.

As shown in FIG. 2, the known OFDM system 10 utilizes an overlapping orthogonal multicarrier modulation technique having a plurality of subcarriers 50. FIG. 3 shows the orthogonal nature of the subcarriers. More particularly, each of four subcarriers 60 of one OFDM data symbol has an integral number of cycles in the interval T. The number of cycles between adjacent subcarriers differs by one.

In one known OFDM transmission system, the complexity of the space-time processor increases as the bandwidth increases and the performance of the space-time processor is significantly degraded when estimated channel parameters are used to construct space-time processors.

It would, therefore, be desirable to provide a MIMO OFDM system having enhanced signal detection. It would further be desirable to increase the accuracy of channel parameter estimates.

SUMMARY OF THE INVENTION

The present invention provides a MIMO OFDM system utilizing independent space-time codes for multiple transmit antenna sets. The receiver decodes the space-time codes using pre-whitening followed by maxiumum likelihood decoding. With this arrangement, the MIMO OFDM system provides spectrally efficient wideband communication. The MIMO-OFDM system determines and uses channel delay profile estimates to achieve more accurate channel parameter estimation.

In one aspect of the invention, a MIMO OFDM communication system includes a plurality of transmit antennas and a plurality of receive antennas. First and second data blocks are each transformed into two signals by respective first and second space-time encoders. Each of the four resultant signals form an OFDM block that is transmitted by a respective transmit antenna. Each receive antenna receives a signal that is the superposition of the four transmitted OFDM blocks. When detecting and decoding the coded signal for the first data block, the coded signal for the second data block is treated as an interfering signal. Similarly, the first data block signal is treated as an interfering signal when detecting and decoding the second data block signal. The received signals are pre-whitened prior to maximum likelihood decoding, which can include Viterbi decoding. In one embodiment, pre-whitening for the maximum likelihood decoding includes minimum mean square error (MMSE) restoration of the desired signals followed by whitening of the residual interfering signals and noise.

Successive interference cancellation can improve performance of the system. More particularly, after the first and second data blocks are decoded, it is determined whether the decoded signals contain errors. In the case where one decoded data block has an error and the other data block does not contain an error, the correct data block signal is regenerated and removed from the received signal. The other data block is then re-detected and decoded from the modified signal.

In a further aspect of the invention, a MIMO OFDM system enhances channel parameter estimation by utilizing relatively accurate channel delay profiles derived from the spatial correlation of the channel impulse responses. In general, the estimated channel responses are weighted based upon a deviation from the average channel response. By more heavily weighting relatively accurate channel responses, the channel parameter estimations are more accurate.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention will be more fully understood from the following detailed description taken in conjunction with the accompanying drawings, in which:

FIG. 1 is a prior art orthogonal frequency division multiplexing (OFDM) system;

FIG. 2 is a prior art graphical depiction showing subchannels used in the OFDM system of FIG. 1;

FIG. 3 is a prior art graphical depiction showing orthogonal subcarriers used in the OFDM system of FIG. 1;

FIG. 4 is a block diagram of a portion of a MIMO-OFDM system in accordance with the present invention;

FIGS. 5A-6B are graphical representations of MIMO-OFDM systems in accordance with the present invention having various channel delay profiles, numbers of receive antennas, and detection techniques;

FIGS. 7A-B are graphical comparisons of MIMO-OFDM systems in accordance with the present invention with ideal and estimated channel parameters; and

FIGS. 8A-D are graphical representations of MIMO-OFDM systems in accordance with the present invention with various Doppler frequencies.

DETAILED DESCRIPTION OF THE INVENTION

In general, the present invention provides an orthogonal frequency division multiplexing (OFDM) system having multiple transmit and receive antennas for forming MIMO channels. With this arrangement, channel estimation and signal detection is enhanced.

FIG. 4 shows a MIMO-OFDM system 100 having multiple, here shown as four, transmit antennas TA1-4 and a plurality of receive antennas RA1-P. Although the MIMO-OFDM system is shown having four transmit antennas, it is understood that any number of transmit antennas can be used. In addition, the number of receive antennas should be equal to or greater than the number of transmit antennas.

The MIMO-OFDM system 100 includes a first space time encoder STE1 that receives a first data block b₁[n,k] and a second space-time encoder STE2 that receives a second data block b₂[n,k]. At time n at tone k, each of the two data blocks, {b_(i)[n,k]: k=0, 1, . . . } for i=1 and 2, is transformed into two signals, {t_(2i+j)[n,k]: k=0, 1, . . . , & j=1, 2} for i=1 and 2, respectively, (Equations 1-3) through the first and second space-time encoders STE1,STE2. Each of the coded signals forms an OFDM block. The transmit antennas TA1-4 transmit the OFDM signals after respective inverse fast Fourier transform IFFT1-4 modulation by respective signals tm_(i)[n,k] for i=1, . . . , 4.

The signals sent by the transmit antennas TA1-4 are received by the receive antennas RA1-RAP. The received signals r₁[n,k], r₂[n,k], . . . , r_(P)[n,k] are transformed by respective fast Fourier transform (FFT) subsystems FFT1-FFTP to generate signals that are provided to a space-time processor STP, which provides detected signal information to respective first and second space-time decoders STD1,STD2. A channel parameter estimator CPE receives the transformed signals from which channel parameter information is determined and then provided to the space-time processor STP for use in decoding the signals.

The received signal at each receive antenna RA1-RAP is the superposition of the four distorted transmitted signals, which can be expressed below in Equation 1:

$\begin{matrix} {{{r_{j}\left\lbrack {n,k} \right\rbrack} = {{\sum\limits_{i = 1}^{4}{{H_{ij}\left\lbrack {n,k} \right\rbrack}{t_{i}\left\lbrack {n,k} \right\rbrack}}} + {w_{j}\left\lbrack {n,k} \right\rbrack}}},} & {{Eq}.\mspace{14mu}(1)} \end{matrix}$ for j=1, . . . , p, where p corresponds to the number of receive antennas, H_(ij)[n,k] denotes the channel frequency response for the k-th tone at time n, corresponding to the i-th transmit and the j-th receive antenna, and w_(j)[n,k] denotes the additive complex Gaussian noise on the j-th receive antenna. The noise is assumed to be zero-mean with variance σ_(n) ² and uncorrelated for different times n's, tones k's, and receive antennas j's.

The input-output relationship for the OFDM signals can be also expressed in vector form as set forth in Equation 2-4 below: r[n, k]=H ₁ [n, k]t ₁ [n, k]+H ₂ [n, k]t ₂ [n, k]+w[n, k],   Eq. (2) where

$\begin{matrix} {{{r\left\lbrack {n,k} \right\rbrack}\overset{\bigtriangleup}{=}\begin{pmatrix} {r_{1}\left\lbrack {n,k} \right\rbrack} \\ \vdots \\ {r_{4}\left\lbrack {n,k} \right\rbrack} \end{pmatrix}},\mspace{14mu}{{w\left\lbrack {n,k} \right\rbrack}\overset{\bigtriangleup}{=}\begin{pmatrix} {w_{1}\left\lbrack {n,k} \right\rbrack} \\ \vdots \\ {w_{4}\left\lbrack {n,k} \right\rbrack} \end{pmatrix}},{{t_{i}\left\lbrack {n,k} \right\rbrack}\overset{\bigtriangleup}{=}\begin{pmatrix} {t_{{2i} + 1}\left\lbrack {n,k} \right\rbrack} \\ {t_{{2i} + 2}\left\lbrack {n,k} \right\rbrack} \end{pmatrix}}} & {{Eq}.\mspace{14mu}(3)} \\ {and} & \; \\ {{H_{i}\left\lbrack {n,k} \right\rbrack}\overset{\bigtriangleup}{=}\begin{pmatrix} {H_{{2i} + 1}{\,_{1}\left\lbrack {n,k} \right\rbrack}} & {H_{{2i} + 2}{\,_{1}\left\lbrack {n,k} \right\rbrack}} \\ \vdots & \vdots \\ {H_{{2i} + 1}{\,_{p}\left\lbrack {n,k} \right\rbrack}} & {H_{{2i} + 2}{\,_{p}\left\lbrack {n,k} \right\rbrack}} \end{pmatrix}} & {{Eq}.\mspace{14mu}(4)} \end{matrix}$

To achieve transmit diversity gain and detection of the transmitted signals, the space-time processor STP extracts the required signals for decoding by the first and second space-time decoders STD1, STD2. The space-time processor and space-time decoders each require channel state information.

In one embodiment, the CPE utilizes conventional training sequences to exploit time and frequency domain correlation of the channel parameters. Exemplary channel estimation techniques are described in Y. Li, et al., “Channel Estimation for OFDM Systems with Transmitter Diversity in Mobile Wireless Channels,” IEEE Journal of Selected Areas in Comm., Vol. 17, pp. 461-471, March 1999, and Y. Li, “Simplified Channel Estimation for OFDM Systems with Multiple Transmit Antennas,” IEEE Journal of Selected Areas in Comm., Vol. 1, Jan., 2002, pp. 67, et seq., which are incorporated herein by reference.

The complex baseband representation of a mobile wireless channel impulse response can be described by Equation 5 below:

$\begin{matrix} {{h\left( {t,\tau} \right)} = {\sum\limits_{k}{{\gamma_{k}(t)}{c\left( {\tau - \tau_{k}} \right)}}}} & {{Eq}.\mspace{14mu}(5)} \end{matrix}$ where τ_(k) is the delay of the k-th path, γ_(k)(t) is the corresponding complex amplitude, and c(t) is a shaping pulse having a frequency response, which can be a square-root raised-cosine Nyquist filter. Due to the motion of the vehicle, the complex amplitudes γ_(k)(t) are wide-sense stationary (WSS), narrow-band complex Gaussian processes, which are independent for each path. The average powers of the complex amplitudes γ_(k)(t) depend on the channel delay profiles, e.g., hilly terrain (HT) and typical urban (TU), which are determined by the environment. The channels typically have the same delay profiles.

From Equation 5 above, it follows that the frequency response at time t can be expressed as set forth below in Equation 6:

$\begin{matrix} {\begin{matrix} {{H\left( {t,f} \right)}\overset{\Delta}{=}{\int_{- \infty}^{+ \infty}{{h\left( {t,\tau} \right)}{\mathbb{e}}^{{- {j2\pi}}\; f\;\tau}{\mathbb{d}\tau}}}} \\ {{= {{C(f)}{\sum\limits_{k}{{\gamma_{k}(t)}{\mathbb{e}}^{{- {j2\pi}}\; f\;\tau_{k}}}}}},} \end{matrix}{where}} & {{Eq}.\mspace{14mu}(6)} \\ {{C(f)}\overset{\Delta}{=}{\int_{- \infty}^{+ \infty}{{c(\tau)}{\mathbb{e}}^{{- {j2\pi}}\; f\;\tau}{\mathbb{d}\tau}}}} & {{Eq}.\mspace{14mu}(7)} \end{matrix}$ In an OFDM system with proper cyclic extension and timing, the channel frequency response, with tolerable leakage, can be represented by Equation 8 below:

$\begin{matrix} {{{H\left\lbrack {n,k} \right\rbrack}\overset{\Delta}{=}{{H\left( {{nT}_{f},{k\;\Delta\; f}} \right)} = {\sum\limits_{l = 0}^{K_{o} - 1}{{h\left\lbrack {n,l} \right\rbrack}W_{K}^{kl}}}}},} & {{Eq}.\mspace{14mu}(8)} \end{matrix}$ where the channel response is

${{h\left\lbrack {n,l} \right\rbrack}\hat{=}{h\left( {{nT}_{f},{k\frac{T_{S}}{K}}} \right)}},{W_{K} = {\exp\left( {{- j}\frac{2\pi}{K}} \right)}},$ K is the number of tones in an OFDM block, T_(f) and Δf are the block length and tone spacing, respectively, and T_(s) is the symbol duration of OFDM, which is related to Δf by T_(s)=1/Δf. The channel responses h[n,1], for l=0, 1, . . . , K₀-1, are WSS, narrow-band complex Gaussian processes. The average power of the channel h[n,l] and index K₀(<K) depend on the delay profiles, e.g., HT and TU, of the wireless channels.

In accordance with the present invention, signal detection is enhanced by utilizing spatial pre-whitening followed by maximum-likelihood (ML) decoding. While joint detection of multiple inputs may be optimal, the concomitant computational complexity renders such an approach impractical.

In contrast to joint detection, in accordance with the present invention the coded signals for the second data block b₂[n,k] are treated as interfering signals when detecting and decoding the first data block b₁[n,k] signals. Similarly, the coded signals for the first data block b₁[n,k] are treated as interfering signals when detecting and decoding the second data block b₂[n,k]. In addition, the transmitted signals are considered to be uncorrelated and Gaussian processes.

For example, the second data block signal b₂[n,k], which is transmitted by the third and fourth transmit antennas TA3,TA4, is treated as an interfering signal when detecting and decoding the first data block b₁[n,k]. Therefore, the interfering signal v[n, k] at the receivers is equal to the channel response vector H₂[n,k] multiplied by the transmitted signal vector t₂[n,k] for the interfering signal plus Gaussian noise w[n,k] as set forth below in Equation 9: v[n, k]=H ₂ [n, k]t ₂ [n, k]+w(n, k)   Eq. (9) The vectors in Equation 9 are derived from Equations 2-4 above.

If v[n, k] is spatially and temporally white, then the minimum Euclidian distance decoder is equivalent to maximum likelihood (ML) decoder. However, if v[n,k] is spatially or temporally correlated, then pre-whitening is used for the ML decoder. Pre-whitening for ML decoding is well known to one of ordinary skill in the art and is described below. In general, the space-time codes for the first and second data blocks are decoded using pre-whitening for flattening the receiver response to facilitate ML detection. The detected constellation points are mapped to the sequence nearest in Euclidean distance.

From the ML criterion, ML decoding is equivalent to finding the estimated first data block {{circumflex over (b)}₁[n,k] that minimizes the value of Equation 10 below:

$\begin{matrix} {{C\left( \left\{ {b_{1}\left\lbrack {n,k} \right\rbrack} \right\} \right)} = {\sum\limits_{k = 1}^{K}{{m\left\lbrack {n,k} \right\rbrack}.}}} & {{Eq}.\mspace{14mu}(10)} \end{matrix}$ where the minimum Euclidean distance is defined in Equations 11-12 below: m[n, k]

(r[n, k]−H₁[n, k]t₁[n, k])^(H)R_(v) ⁻¹[n, k](r[n, k]−H₁[n, k]t₁[n, k])   Eq. (11) and R_(v)[n, k]

E{v[n, k]v^(H)[n, k]}=H₂[n, k]H₂ ^(H)[n, k]+σ_(n) ²I   Eq. (12) As known to one of ordinary skill in the art, the inverse of the matrix R_(v)[n,k], which is derived from the cross correlation H₂[n,k]H₂ ^(H)[n,k] of the interfering channel response, can be used to whiten the signal.

Direct calculation then yields that the minimum Euclidean distance can be defined as set forth in Equation 13

$\begin{matrix} {{m\left\lbrack {n,k} \right\rbrack} = {{{r^{H}\left\lbrack {n,k} \right\rbrack}{R_{v}^{- 1}\left\lbrack {n,k} \right\rbrack}{r\left\lbrack {n,k} \right\rbrack}} - {{t_{1}^{H}\left\lbrack {n,k} \right\rbrack}{H_{1}^{H}\left\lbrack {n,k} \right\rbrack}{R_{v}^{- 1}\left\lbrack {n,k} \right\rbrack}{r\left\lbrack {n,k} \right\rbrack}} - {{r^{H}\left\lbrack {n,k} \right\rbrack}{R_{v}^{- 1}\left\lbrack {n,k} \right\rbrack}{H_{1}\left\lbrack {n,k} \right\rbrack}{t_{1}\left\lbrack {n,k} \right\rbrack}} + {{t_{1}^{H}\left\lbrack {n,k} \right\rbrack}{H_{1}^{H}\left\lbrack {n,k} \right\rbrack}{R_{v}^{- 1}\left\lbrack {n,k} \right\rbrack}{H_{1}\left\lbrack {n,k} \right\rbrack}{{t_{1}\left\lbrack {n,k} \right\rbrack}.}}}} & {{Eq}.\mspace{14mu}(13)} \end{matrix}$ Denoting {tilde over (H)}[n,k] as a two by two matrix satisfying Equation 14 and 15 H₁ ^(H)[n, k]R_(v) ⁻¹[n, k]

{tilde over (H)}^(H)[n, k]{tilde over (H)}[n, k]  Eq. (14) and L₁[n, k]

(H₁{tilde over (H)}₁ ⁻¹)^(H)R_(v) ⁻¹[n, k].   Eq. (15) it then follows that the minimum Euclidean distance can be represented in Equation 16

$\begin{matrix} \begin{matrix} {{m\left\lbrack {n,k} \right\rbrack} = {{{r^{H}\left\lbrack {n,k} \right\rbrack}{R_{v}^{- 1}\left\lbrack {n,k} \right\rbrack}{r\left\lbrack {n,k} \right\rbrack}} -}} \\ {{{t_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}^{H}\left\lbrack {n,k} \right\rbrack}{L_{1}\left\lbrack {n,k} \right\rbrack}{r\left\lbrack {n,k} \right\rbrack}} -} \\ {{{r^{H}\left\lbrack {n,k} \right\rbrack}{L_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}\left\lbrack {n,k} \right\rbrack}{t_{1}\left\lbrack {n,k} \right\rbrack}} +} \\ {{t_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}\left\lbrack {n,k} \right\rbrack}{t_{1}\left\lbrack {n,k} \right\rbrack}} \\ {= {{{r^{H}\left\lbrack {n,k} \right\rbrack}{R_{v}^{- 1}\left\lbrack {n,k} \right\rbrack}{r\left\lbrack {n,k} \right\rbrack}} -}} \\ {{{t_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{r}}_{1}\left\lbrack {n,k} \right\rbrack}} -} \\ {{{{\overset{\_}{r}}_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}\left\lbrack {n,k} \right\rbrack}{t_{1}\left\lbrack {n,k} \right\rbrack}} +} \\ {{t_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{H}}_{1}\left\lbrack {n,k} \right\rbrack}{t_{1}\left\lbrack {n,k} \right\rbrack}} \\ {= {{{r^{H}\left\lbrack {n,k} \right\rbrack}{R_{v}^{- 1}\left\lbrack {n,k} \right\rbrack}{r\left\lbrack {n,k} \right\rbrack}} - {{{\overset{\_}{r}}_{1}^{H}\left\lbrack {n,k} \right\rbrack}{{\overset{\_}{r}}_{1}\left\lbrack {n,k} \right\rbrack}} +}} \\ {{{{{\overset{\_}{r}}_{1}\left\lbrack {n,k} \right\rbrack} - {{{\overset{\_}{H}}_{1}\left\lbrack {n,k} \right\rbrack}{t_{1}\left\lbrack {n,k} \right\rbrack}}}}^{2},} \end{matrix} & {{Eq}.\mspace{14mu} 16} \end{matrix}$ where {tilde over (r)}₁[n, k]

L₁[n, k]r[n, k].   Eq. (17)

When the well known Viterbi algorithm is used for the ML decoding of the space-time codes, the first two terms in the above equation for m[n,k] are independent of the detected data and only the third term, i.e., ∥{tilde over (r)}₁[n, k]−{tilde over (H)}₁[n, k]t₁[n, k]∥² , is related to the detected data and affects the metric in the trellis search when the Viterbi algorithm is used. Trellis searching is well known to one of ordinary skill in the art. Therefore, the ML decoding is equivalent to finding the estimated first data block {{circumflex over (b)}₁[n, k] that minimizes the value of Equation 18 below:

$\begin{matrix} {{\overset{\_}{C}\left( \left\{ {b_{1}\left\lbrack {n,k} \right\rbrack} \right\} \right)} = {\sum\limits_{k = 1}^{K}{{{{\overset{\_}{r}}_{1}\left\lbrack {n,k} \right\rbrack} - {{{\overset{\_}{H}}_{1}\left\lbrack {n,k} \right\rbrack}{t_{1}\left\lbrack {n,k} \right\rbrack}}}}^{2}}} & {{Eq}.\mspace{14mu}(18)} \end{matrix}$ Thus, after pre-whitening, a conventional space-time decoder for a 2-transmit and 2-receive antenna system can be used.

Note that L₁[n,k] can be rewritten as shown in Equation 19: L ₁ [n, k]=({tilde over (H)} ⁻¹ [n, k])^(H) H ₁ ^(H) [n, k]R _(v) ⁻¹ [n, k].   Eq. (19) A predetermined weight matrix H₁ ^(H)[n, k]R_(v) ⁻¹[n, k] for minimum mean-square error (MMSE) restoration of the desired signal t₁[n,k] suppresses the interfering signal t₂[n,k]. After MMSE signal restoration, the correlation matrix of the residual interferers and noise can be represented as set forth in Equation 20: E{H ₁ ^(H) [n, k]R _(v) ⁻¹ [n, k]v[n, k](H ₁ ^(H) [n, k]R _(v) ⁻¹ [n, k]v[n, k])^(H) }=H ₁ ^(H) [n, k]R _(v) ⁻¹ [n, k]H ₁ [n, k]. The term {tilde over (H)}⁻¹[n, k]^(H) whitens the residual interferers and noise. Therefore, the pre-whitening processing for the ML decoder includes MMSE restoration of the desired signals followed by whitening of the residual interferers and noise.

In another aspect of the invention, successive interference cancellation is utilized to improve the overall system performance. The successive interference cancellation can be based upon cyclic redundancy check (CRC) codes and/or signal quality.

In one embodiment, a MIMO-OFDM system includes successive interference cancellation based upon CRC codes for identifying decoding errors. First and second data blocks b₁[n,k], b₂[n,k] can be decoded as described above. If an error is detected in one of the data blocks, such as the first data block b₁[n,k], and no errors are detected in the other data block, such as the second data block b₂[n,k], then the coded signals for the correct (second) data block are regenerated at the receiver and removed from the received signal. That is, the coded signals for second data block b₂[n,k] are removed from the received signal such that cleaner signals can be used to re-detect and decode the first data block b₁[n,k], which originally contained an error. More particularly, the first data block can be re-detected and decoded without interference from the second data block signals.

In a further embodiment, a MIMO-OFDM system detects and decodes each of the data block signals, which do not include error codes such as CRC codes. As is known to one of ordinary skill in the art, some systems do not utilize error coding when signal quality is expected to be above a predetermined threshold. In accordance with the present invention, each decoded signal has an associated MMSE. The signal with the higher quality, e.g., lower MMSE, is removed such that it does not interference with the other signal.

In another aspect of the invention, a MIMO-OFDM system includes enhanced channel parameter estimation using relatively accurate channel delay profiles. In one embodiment, a known decision-directed channel parameter estimator along with optimum training sequences for OFDM systems with multiple transmit antennas are used to estimate the channel delay profiles. Illustrative parameter estimations and training sequences are described in Y. (Geoffrey) Li, et al. “Channel Estimation for OFDM Systems with Transmitter Diversity and its Impact on High-Rate Data Wireless Networks,” IEEE Journal of Selected Areas in Comm., Vol. 17, pp. 461-471, March 1999, and Y. Li, et al., “Simplified Channel Estimation for OFDM Systems with Multiple Transmit Antennas,”mentioned above. The estimated channel delay profiles are used to determine channel parameter estimations.

In general, the time and frequency correlation of the channel profiles are used to more accurately estimate the channel parameter estimation. More particularly, channel profiles are correlated in time since the change over time is relatively slow. Similarly, adjacent channels have similar frequency responses and are thus correlated with respect to frequency. By determining the average channel impulse response for example, a deviation from the average can form the basis to weight each channel impulse response. With this weighting arrangement, the channel delay estimates can be more accurately estimated to enhance channel parameter estimation.

The channel impulse response h_(ij)[n,l] can be estimated using the correlation of channel parameters in the time and frequency domains. The estimated channel impulse response ĥ_(ij)[n,l] can be reconstructed using a discrete Fourier transform (DFT) as shown below in Equation 21

$\begin{matrix} {{{\hat{H}}_{ij}\left\lbrack {n,k} \right\rbrack} = {\sum\limits_{l = 0}^{K_{o} - 1}{{{\hat{h}}_{ij}\left\lbrack {n,l} \right\rbrack}W_{K}^{kl}}}} & {{Eq}.\mspace{14mu}(21)} \end{matrix}$ where ĥ_(ij)[n,l] contains the true channel parameter h_(ij)[n,l], W_(k) ^(kl) is described above in Equation 5, and an estimation error e_(ij)[n,l] as expressed below in Equation 22: ĥ _(ij) [n, l]=h _(ij) [n, l]+e _(ij) [n, l]  Eq. (22) The estimation error e_(ij)[n,l] can be Gaussian with zero-mean and variance σ². The parameter estimation quality can be measured with a normalized MSE (NMSE) defined as follows in Equation 23:

$\begin{matrix} {{N\; M\; S\; E}\overset{\Delta}{=}\frac{E{{{{\hat{H}}_{ij}\left\lbrack {n,k} \right\rbrack} - {H_{ij}\left\lbrack {n,k} \right\rbrack}}}^{2}}{E{{H_{ij}\left\lbrack {n,k} \right\rbrack}}^{2}}} & {{Eq}.\mspace{14mu}(23)} \end{matrix}$ Thus, the NMSE for the channel response can be expressed as set forth in Equation 24: NMSE_(r)=K_(o)σ²   Eq. (24) assuming that the NMSE is normalized as shown in Equation 25:

$\begin{matrix} {{\sum\limits_{l = 0}^{K_{o} - 1}{E{{h_{ij}\left\lbrack {n,l} \right\rbrack}}^{2}}} = {{\sum\limits_{l = 0}^{K_{o} - 1}\sigma_{l}^{2}} = 1}} & {{Eq}.\mspace{14mu}(25)} \end{matrix}$ with σ_(l) ²

E∥h_(ij)[n, l]∥².

If the channel delay profile is known, that is, σ₁ ² for l=0, . . . , K₀-1 is known, and is used to reconstruct channel frequency response from ĥ_(ij)[n, l], the MSE of Ĥ_(ij)[n, k], can be significantly reduced. In this case, if weighting factors α₁'s are selected to minimize the NMSE of the estimated channel response as in Equation 26:

$\begin{matrix} {{{\hat{H}}_{ij}\left\lbrack {n,k} \right\rbrack}\overset{\Delta}{=}{\sum\limits_{l = 0}^{K_{o} - 1}{\alpha_{l}{{\hat{h}}_{ij}\left\lbrack {n,l} \right\rbrack}W_{K}^{kl}}}} & {{Eq}.\mspace{14mu}(26)} \end{matrix}$ then the optimal weighting factor α₁ can be defined in Equation 27:

$\begin{matrix} {\alpha_{l} = \frac{\frac{\sigma_{l}^{2}}{\sigma_{l}^{2} + \sigma^{2}}}{\sum\limits_{m = 0}^{K_{o} - 1}\frac{\sigma_{m}^{4}}{\sigma_{m}^{2} + \sigma^{2}}}} & {{Eq}.\mspace{14mu}(27)} \end{matrix}$

which can also be expressed as an inverse of the sum of one plus a ratio of noise power to a power of channel response estimates averaged over the receive antennas, i.e.,

$\frac{1}{1 + \frac{\sigma^{2}}{\sigma_{l}^{2}}}$

and the resultant NMSE is expressed in Equation 28:

$\begin{matrix} {{NMSE}_{o} = \frac{\sigma^{2}{\sum\limits_{m = 0}^{K_{o} - 1}\frac{\sigma_{m}^{2}}{\sigma_{m}^{2} + \sigma^{2}}}}{\sum\limits_{m = 0}^{K_{o} - 1}\frac{\sigma_{m}^{4}}{\sigma_{m}^{2} + \sigma^{2}}}} & {{Eq}.\mspace{14mu}(28)} \end{matrix}$

Channel delay profiles depend upon the environment and therefore are usually unknown to users of the system. However, for MIMO-OFDM systems, channels corresponding to different transmit or receive antennas should have the same delay profiles. Thus, σ₁ ²=E∥h_(ij)[n , l]∥² can be estimated by averaging the estimated channel impulse response for the channels formed by the four transmit antennas and the p receive antennas as shown below in Equation 29:

$\begin{matrix} {{\hat{\sigma}}_{l}^{2} = {\frac{1}{4p}{\sum\limits_{i = 1}^{4}{\sum\limits_{j = 1}^{p}{{{{\hat{h}}_{ij}\left\lbrack {n,l} \right\rbrack}}^{2}.}}}}} & {{Eq}.\mspace{14mu}(29)} \end{matrix}$

With the estimated variance, channel delay profile estimates are more accurate so as to improve channel parameter estimation.

EXAMPLE

An exemplary MIMO-OFDM system in accordance with the present invention was simulated. The known typical urban (TU) and hilly terrain (HT) delay profiles were used with Doppler frequencies of 5, 40, 100, and 200 Hz, respectively. The channels corresponding to different transmit or receive antennas have the same statistics. Four transmit antennas and different numbers of receive antennas were used to form a 4-input multiple-output OFDM systems.

To construct an OFDM signal, the entire channel bandwidth, 1.25 MHz, is divided into 256 subchannels. The 2 subchannels on each end are used as guard tones, and the remaining (252 tones) are used to transmit data. To make the tones orthogonal to each other, the symbol duration is about 204.8 microseconds. An additional 20.2 microseconds guard interval is used to provide protection from intersymbol interference due to channel multipath delay spread. This results in a total block length T_(f) of about 225 microseconds and a subchannel symbol rate r_(b) of 4.44 kbaud.

A 16-state space-time code with 4-PSK is used. Each data block, containing 500 bits, is coded into two different blocks, each of which has exactly 252 symbols, to form an OFDM block. Therefore, the OFDM system with 4 transmit antennas can transmit 2 space-time codewords (1000 bits in total). Each time slot includes 10 OFDM blocks, with the first block used for training and the remaining nine blocks used for data transmission. Thus, the system can transmit data at a rate of 4 Mbits/sec over a 1.25 MHz channel, i.e., the transmission efficiency is 3.2 bits/sec/Hz.

FIGS. 5A-6B show the simulated performance of MIMO-OFDM systems in accordance with the present invention with different channel delay profiles, numbers of receive antennas, and detection techniques. FIG. 5A demonstrates the effect of interleaving on the performance improvement. By interleaving, the required SNR for a 10% WER is improved by 1.5 dB for the TU channel and by 0.7 dB for the HT channel. Since the HT channel has more diversity than the TU channel before interleaving, interleaving has less gain for the HT channel than for the TU channel.

FIG. 6A-B compare the WERs of a system with interleaving for different detection techniques. As shown, the system with successive interference cancellation based on CRC and signal quality (MMSE) can reduce the required WERs for a 10% WER by 2.5 and 1.8 dB, respectively. All the performance curves in FIGS. 5A-5B are for OFDM with 4 transmit and 4 receive antennas. As the number of receive antennas increases, performance improves, as can be seen from FIG. 6A-6B. In particular, if the receive antenna number is increased from 4 to 6, the OFDM system with either the TU or HT channels will have about 4 dB performance improvement.

FIGS. 7A-7B compare the performance of MIMO-OFDM systems with the ideal or the estimated channel parameters for different channels with a 40 Hz Doppler frequency. From FIG. 7A it can be seen that the required SNRs for a 10% WER are 10-11 dB for a MIMO-OFDM system with estimated channel parameters for successive interference suppression and space-time decoding, which is 1.5-2 dB higher than with ideal channel parameters for signal detection and decoding. With additional receive antennas, performance improves, as shown in FIG. 7B. In particular, for a system with estimated channel parameters, the required SNR for a 10% WER is reduced by 4.5 dB and 2 dB when the number of receive antennas is increased from 4 to 6 and 6 to 8, respectively.

FIGS. 8A-D compare the performance of OFDM systems with different Doppler frequencies. As the Doppler frequency becomes higher, the channel estimation error increases, and therefore the system suffers more degradation. For a MIMO-OFDM system with 4 transmit and 4 receive antennas, the required SNR for a 10% WER is degraded by 2.4 dB when the Doppler frequency is increased from 40 Hz to 100 Hz. However, with more receive antenna numbers, the degradation is reduced. The degradation is only about 0.4 dB with 10 receive antennas.

The present invention provides an OFDM system having multiple transmit and receive antennas to form a multiple-input multiple-output (MIMO) system that increases system capacity. A pre-whitening technique for ML decoding and successive interference cancellation technique are disclosed. Using these techniques in a four input/four output OFDM system, the net data transmission rate can reach 4 Mbits/sec over a 1.25 MHz wireless channel with a 10-12 dB SNR required for a 10% WER, depending on the radio environment and signal detection technique for word lengths up to 500 bits. A MIMO-OFDM system in accordance with the present invention can be effectively used in high data-rate wireless systems

One skilled in the art will appreciate further features and advantages of the invention based on the above-described embodiments. Accordingly, the invention is not to be limited by what has been particularly shown and described, except as indicated by the appended claims. All publications and references cited herein are expressly incorporated herein by reference in their entirety. 

What is claimed is:
 1. A receiver constructed to receive signals from a remote unit that at time n transmits a data block b_(i)[n,k] over two transmit antennas and transmits a data block b₂[n,k] over two other transmit antennas, where k is a tone and k=0, 1, 2, . . . , comprising: N receive antennas, each receiving a signal that includes signals related to the signals transmitted by the transmit antennas of the remote unit, where N is equal or greater than 4; N fast Fourier transform (FFT) modules, each coupled to an associated one of said receive antennas; a space-time processing (STP) unit, responsive to said FFT modules, to pre-whiten signals supplied to said STP unit and to perform maximum likelihood detection on each of the pre-whitened signals, and which, when initially detecting data pertaining to block b₁[n,k] treats other received signals, including signals received from the antennas that transmit data pertaining to block b₂[n,k], as interfering noise, and when initially detecting data pertaining to block b₂[n,k] treats other received signals, including signals received from the antennas that transmit data pertaining to block b₁[n,k], as interfering noise; a first space-time decoder responsive to output signals of said STP unit for decoding block b₁[n,k] ; a second space-time decoder responsive to output signals of said STP unit for decoding block b₂[n,k]; and a channel parameter estimator responsive to output signals of said FFT modules and output signals of said space-time decoders, for developing parameter estimates of channels between the remote unit and the receiver, and applying the developed parameter estimates to said space-time decoders.
 2. The receiver of claim 1 where estimates of channel delay profiles are used to determine said parameter estimates.
 3. The receiver of claim 1 further comprising a feedback signal path from at least one of the space-time decoders to said STP unit, which enables said STP unit to modify results of said maximum likelihood detection based on signals supplied by said feedback signal path.
 4. The receiver of claim 1 further comprising: a feedback signal path from at least one of the space-time decoders to said STP unit; and when (a) initial detecting by said STP unit of data pertaining to block b₁[n,k] and data pertaining to block b₂[n,k] is completed, (b) said first and second space time decoders responsively form decoded data blocks b₁[n,k] and b₂[n,k], respectively, and (c) if it is determined that an error exists in decoded data block b₁[n,k] and no error is detected in decoded data block b₂[n,k], said STP unit (i) encodes the decoded data block b₂[n,k], which is received via said feedback path, (ii) removes the encoded data block from the signals received by said STP unit from said FFT modules, to form modified signals, and (iii) responsive to the modified signals, detects data pertaining to block b₁[n,k] and data pertaining to block b₂[n,k].
 5. A method comprising: at each of N receive antennas, receiving a signal that includes signals related to signals transmitted from a remote unit at time n, the signals including a data block b₁[n,k] transmitted over first transmit antennas and a data block b₂[n,k] transmitted over second transmit antennas, where k is a tone and k=0, 1, 2, . . . , where N is equal or greater than 4; transforming the received signal in N fast Fourier transform (FFT) modules coupled to an associated one of said receive antennas and generating respective transformed signals; receiving the transformed signals at a space-time processing unit and pre-whitening signals in the space-time processing unit associated with the data blocks b₁[n,k] and b₂[n,k] and performing maximum likelihood detection on each of the pre-whitened signals, and when initially detecting data pertaining to block b₁[n,k] treating other received signals, including signals received from the second transmit antennas that transmit data pertaining to block b₂[n,k], as interfering noise, and when initially detecting data pertaining to block b₂[n,k] treating other received signals, including signals received from the first transmit antennas that transmit data pertaining to block b₁[n,k], as interfering noise; decoding block b₁[n,k] in a first space-time decoder; decoding block b₂[n,k] in a second space-time decoder; developing channel parameter estimates using output signals of said FFT modules and output signals of said space-time decoders; and using the channel parameter estimates in said first and second space-time decoders.
 6. The method of claim 5 further comprising using estimates of channel delay profiles to determine said channel parameter estimates.
 7. The method of claim 5 further comprising modifying results of said maximum likelihood detection based on signals from at least one of the space-time decoders.
 8. The method of claim 5 further comprising: said first and second space time decoders forming decoded data blocks b₁[n,k] and b₂[n,k], respectively, responsive to completion of initial detecting by said STP unit of data pertaining to block b₁[n,k] and data pertaining to block b₂[n,k]; if it is determined that an error exists in decoded data block b₁[n,k] and no error is detected in decoded data block b₂[n,k], said STP unit encoding the decoded data block b₂[n,k], which is received via a feedback path from at least one of the space-time decoders to said STP unit, removing the encoded data block from the signals received by said STP unit from said FFT modules, to form modified signals, and responsive to the modified signals, detecting data pertaining to block b₁[n,k]. 